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I. INTRODUCTION 

The A(3872) was discovered at Belle [l[ and then later 
was also observed at CDFII and DO collaborations and 
BaBar 0, In all these experiments the X has been 

discovered and observed in the decay channel J/yjir + n~ . 
There is strong evidence that the dipion generated in this 
decay channel comes from a p meson |5| . Later on also the 
decays of the X into J / 'ipir + ir~ w° and J/ip"f have been 
observed @ , this latter decay channel indicating that the 
C-parity of the X is positive. The quantum numbers of 
the A" (3872) have been investigated in [7(, concluding 
that it must correspond to J p = 1 ++ or J p = 2 h . Ob- 
served only on a neutral charge state it is assumed to have 
isospin 1 = 0. Its decay into J /ybr] has been investigated 
in [8j but only an upper bound has been found. The non 
observation of the decay J/tpr] is a further evidence of the 
positive C-parity of the X. It could also mean that in 
the particular reaction of [8[ the A (3872) was necessar- 
ily produced with positive C-parity, without ruling out 
the possibility of a nearby state with negative C-parity. 
The existence of two nearly degenerate A (3872) states 
appears in some theoretical models [H, E3|- The most 
popular view about the nature of this resonance is that 
it is made of DD* Eljp, H El EI El El , a recent 
review can be seen in 17J • One of the problems faced by 
these models is the large ratio for 



B(X -> J/yjir+n-n ) 
BjX J/t/m+n-) 



= 1.0 ±0.4 ±0.3 



(1) 



Indeed, since the resonance has positive C-parity the de- 
nominator can go via J/iftp as supported by the experi- 
ment |5(. However, the A(3872) state has 1 = and then 
isospin is violated. On the contrary the numerator can go 
through J/ipLu as supported by experiment d, in which 
case there is no violation of isospin. The fact that the 
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ratio is so large in spite of the violation of isospin found 
a plausible explanation in [11 EB| , where the state was 
supposed to be largely D°D*° but with some coupling to 
both J/ipui and J/ipp. Even if the coupling to J/ipP is 
small, as expected from isospin symmetry breaking, the 
larger phase space for J/tpp decay than for J / tpu, because 
of the large width of the p, can account for the large ratio. 
Although other charged DD* components can appear in 
the wave function, the neutral charge component is pre- 
ferred since it is the one closest to threshold and hence 
should have the largest weight. The idea is intuitive and 
widely accepted, see (T^. The idea on the dominance 
of the neutral component is worth pursuing. Indeed. In 
[l(| , where a dynamical theory for the generation of the 
A (3872) resonance based on the hidden gauge approach 
for the vector-meson interaction was done, isospin sym- 
metry was kept and the masses of the charged and neu- 
tral D mesons were taken equal. The fact that the bind- 
ing energy for the D°D*° is so small advises to revise 
the model to account for the mass differences with the 
charged DD*, which can induce isospin breaking and a 
dominance of the D°D*° in the wave function. In the 
present paper we will face this problem and will discuss 
qualitatively as well as quantitatively, the limit of zero 
binding energy, with interesting results. We will also re- 
visit the interpretation of the reaction production of DD* 
[24| showing that it gives support to the existence of the 
A (3872) as a DD* narrow state of positive parity, with- 
out ruling out the possible existence of a broader one of 
negative parity. We discuss the decay modes of the neg- 
ative C-parity state and speculate on where could it be 
found and the possible difficulties in its observation. 

The present work should also be looked at with some 
perspective. Although most of the work for the A (3872) 
has concentrated on the molecular DD* picture, as we 
have mentioned, since the mass of is so close to the 
threshold, there are other pictures proposed to describe 
it in terms of quarks, tetraquarks and other dynamical 
pictures (see [20, HH, [22j] for reviews). The molecular 
picture of some meson resonances is catching up, par- 
ticularly when it comes to interpret states that do not 
fit clearly in a qq picture. This is the case of the A, Y 
and Z resonances recently discovered at the i3-factories 
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for which there are also other structures proposed and 
about which there is an intensive debate (sec 23] for a 
recent overview on the subject). 

This work is organized as follows: in the next sec- 
tion we shortly explain our phenomenological Lagrangian 
and present our framework to generate dynamically reso- 
nances from the interaction of pseudoscalars with vector- 
mesons. In the same section we show results with and 
without isospin violation and we also show what happens 
when the binding energy goes to zero. In section III wc 
comment on the two C-parity states that our model gen- 
erates and in section IV we make our final remarks and 
conclusions. 



<I> 



Note that these fields differ from those used in [1(3] be- 
cause of the inclusion of r\-rj and u-<f> mixing. 
For each one of these fields a current is defined: 

J M = (d^)t> - <5>d^ (4) 
J„ = id li V v )V-V u d M V v . (5) 

The Lagrangian is constructed by coupling these cur- 
rents: 

£ppw = -^Tr(J^). (6) 

In the way it is constructed this Lagrangian is SU (4) 
symmetric, but we know that SU (4) symmetry is badly 
broken in nature. To take this into account we will 
break the SU(4) symmetry of the Lagrangian in the 
following way: Assuming vector-meson dominance we 
recognize that the interaction behind our Lagrangian 
is the exchange of a vector meson in between the two 
hadronic currents. If the initial and final pseudoscalars 



II. THEORETICAL MODEL 



The framework that we describe here is explained in 
more details in [l(| and references therein. 

We want to study the interaction of a pseudoscalar 
with a vector-meson that, in s-wave, has the quantum 
numbers of an axial: 1 + . The starting point of our model 
consists of the fields belonging to the I5-plet and a singlet 
of SU(4) describing the pseudoscalar and vector-mesons: 



(2) 



(3) 



(and vector-mesons), in a given process, have different 
charm quantum number, it means that the vector-meson 
exchanged in such a process is a charmed meson, and 
hence a heavy one. In these cases we suppress the term 
in the Lagrangian containing such processes by a factor 
7 = m 2 L /m 2 H where is the typical value of a light 
vector-meson mass (800 MeV) and rnu the typical value 
of the heavy vector- meson mass (2050 MeV). We also 
suppress, in the interaction of D-mesons the amount of 
the interaction which is driven by a J /if) exchange by 
the factor rp = m^/m^. Another source of symmetry 
breaking will be the meson decay constant / appearing 
in the Lagrangian. For light mesons we use / = f w = 93 
MeV but for heavy ones / = fr> = 165 MeV. 

So, for a given process (P(p)V(fc))i -► {P' {p')V {k')) 3 
we have the amplitude: 

M tJ (s,t,u) = -_|L-( s -u) e .e' (7) 

where s and u are the usual Mandelstam variables, /, 
is the pseudoscalar i meson decay constant, e are the 
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vector-meson polarization vectors and i, j refer to the 
initial and final channels in the coupled channel space. 
The coefficient matrices £y can be directed calculated 
from the Lagrangian of eq. © in charge basis. The £y 
coefficients in charge basis are given in the appendix. 

The amplitude in eq. JJJ) is projected in s-wave and 
plugged into the scattering equation for the coupled chan- 
nels: 



A. Isospin Symmetric Case 

If we set the masses of all mesons belonging to a same 
isospin multiplet to a common value, our results will be 
isospin symmetric. Moreover we can consider the trans- 
formation under C-parity of pseudoscalar and vector- 
mesons in order to construct C-parity symmetric states: 



T 



V + VGT. 



(8) 



In this equation G is a diagonal matrix with each one of 
its elements given by the loop function for each channel 
in the coupled channel space. For channel i with mesons 
of masses mi and m-i Gu is given by: 



Ga — 



16tt 2 



Log— 



Log- 



2s 



'Log j 



ml + 2py/s 



■2py/aJ 



(9) 



where p is the three momentum of the two mesons in 
the center of mass frame. The two parameters p and a 
are not independent, we fix p=1500 MeV and change a 
to fit our results within reasonable values in the natural 
range [25]. We actually use two a as free parameters, 
one for loops with only light mesons in the channel i. we 
call cxl and set it to a.£,=-0.8, the value used in [10| to 
fit the low lying axial resonances. For loops with heavy 
particles we use ajj and we comment, in what follows, 
the effects of changing this parameter. Note that just 
the combination a — Log(p 2 ) is the free parameter in 
G of eq. (JOj) . The parameter p is there to set a scale 
of energies, it can be chosen arbitrarily and hence is not 
related with the cut off in three momentum which can be 
alternatively used to evaluate Gu from the loop function 
of two meson propagators. Once p is fixed then there 
is a relationship between a of eq. (O and an equivalent 
cut off |26J. In any case, the equivalence of the cut off 
method and the dimensional regularization of eq. ([9]) 
in a certain region of energies requires the cut off to be 
reasonably bigger than the on shell three momenta of the 
intermediate states, which is fulfilled in the energy regime 
that we study here. In the present case, the DD* states, 
with 20 MeV above threshold have a three momentum 
around 195 MeV, while the cut offs are of the order of 
650-850 MeV, as we shall see later on. 

The imaginary part of the loop function ensures that 
the T-matrix is unitary, and since this imaginary part is 
known, it is possible to do an analytic continuation for go- 
ing from the first Riemann sheet to the second one. Pos- 
sible physical states (resonances) arc identified as poles 
in the T-matrix calculated in the second Riemann sheet 
for the channels which have the threshold below the res- 
onance mass. 



CP = P 
CV = -V 



(10) 

(11) 



States like DD* and DD* mix up to form a positive 
and a negative C-parity state. The same happens for the 
kaons and the D s mesons. If one writes the £jj coefficients 
that appear in the amplitude of eq. |(7j) in C-parity basis 
for all two meson states with quantum numbers C=0 and 
S=0, the coupled channel space splits into two, a positive 
and a negative C-parity part. These coefficients are given 
in the appendix. 

In charge basis, for positive C-parity one has the fol- 
lowing channels: K*°K° — c.c, p + ir~ — cc, D*°D° — cc, 
D*+D- - c.c, D*+D- - c.c. and K*+ K~ - c.c. While 
for negative C-parity the channels are: /9 + 7r~ + c.c, 
K*+K- + c.c, p°ir°, W7T°, <f>ir°, p r], p°r]t, K*°K° + c.c, 

D*+D- + CC, D*°D° + C.C, p°T] c , J/-07T , LOT], <f>7), UT]t, 

cjrr)', Lur) c , (pr] c , J/ipr), J/ipr/t, D* + D~ + c.c. and J/ipr] c . 

For the masses of the mesons we use the following val- 
ues: 

771^=137.5 MeV, m K =496 MeV, m v =548 MeV, 
m D =1867.5 MeV, m Ds =1968 MeV, to,, c =2980 MeV, 
to,,,=958 MeV, m p =775 MeV, m K .=894 MeV, 771^=783 
MeV, m =lO19 MeV, m D , =2008.5 MeV, to d -=2II2 
MeV and m J/v ,=3097 MeV. 

If we set ct#=-1.34, which is equivalent to a cut-off of 
830 MeV in the three momentum, we get two poles with 
opposite C-parity, the positive one at 3866 MeV with a 
width smaller than I MeV and the negative one at (3875- 
25i) MeV, which means a width around 50 MeV. The 
poles appear in isospin 1=0, as we determine from com- 
bining the charge states into definite isospin states. Now 
while increasing the value of an (lowering the cut-off) the 
poles approach the threshold (at 3876 MeV in the isospin 
symmetric case) . The negative C-parity pole touches the 
threshold for an values bigger that -1.33 (cut-off of 820 
MeV) , while the positive C-parity one reaches the thresh- 
old for an around -1.185 (cut-off equivalent to 660 MeV). 
Once the pole crosses the threshold it docs not appear in 
the second Riemann sheet, it is no longer a resonance, but 
becomes a virtual state. Yet a peak can be seen in the 
cross section of some channels, but can not be identified 
as a pole in the second Riemann sheet of the T-matrix. 



B. Isospin Breaking 

In order to investigate the isospin breaking we define 
the following quantities: 
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Am T =2.5 MeV, Am K =-2 MeV, Am D =2.5 MeV, 
Am K *=-2 MeV and Am D .=1.5 MeV. In this way the 
masses of the members of a multiplet split: for the 
charged members of a multiplet the mass will be equal to 
771+ Am while for the neutral members it will be m — Am. 

Now there are two DD* thresholds nearby, the neutral 
one at 3872 MeV and the charged one at 3880 MeV. The 
AT (3872) state is a very weakly D°D*° bound state and 
the fact that the binding energy is much smaller than 
the difference between these two thresholds could reflect 
itself in a large isospin violation in observables. 

For simplicity let us consider, for the moment, a toy 
model with only two channels, with neutral and charged 
D and D* mesons. In this model we assume the potential 
V to be a 2x2 matrix: 
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(12) 



with v constant, which indeed is very close to the real 
one in a small range of energies. 

In this case the solution of the scattering equation © 

is: 



T 



V 



1-vGu-vG 



22 



(13) 



where Gn and G22 are the loop function calculated for 
channels 1 and 2 respectively. If there is a pole at s=sr 
we can expand T close to this pole as: 



T ■ — 



sr 



(14) 



where gi is the coupling of the pole to the channel i. 
The product gigj is the residue at the pole and can be 
calculated with: 



to their building blocks in studies [13, [23] of dynamically 
generated states following the method of the composit- 
ness condition of Weinberg [2^, [3(| • 

We will come back to this issue again with the realistic 
model. Now, in what follows, the arguments used do not 
require the toy model any longer. 

Suppose that the V(3872) decays through the diagram 
in figure [2] In this figure the D mesons can be either 
charged or neutral. For the isospin 1=1 state with the 
p meson in the final state, the diagrams with neutral D 
mesons interfere destructively with those with charged D 
mesons, while in the lo case they sum up. If the vertices 
have the same strength for p and u> production (this is 
the case in the framework of the hidden gauge formalism 
[3lL [32l . l33l |) the ratio of the amplitudes will be given 
by the ratio of the difference between the charged and 
neutral loops divided by the sum of the loops: 



R 



p/u 



Gn — G22 
Gn + G22 



(17) 



lim (s - s R )Tij = lim (s - s R )- — ^ 77-f 15 ) 

s^s R 8—>8r 1 — vGn — V(j22 

We can apply the l'Hopital rule to this expression and 
we get: 



lim (s 

S^Sr 



ViA 



<lG- 2 



(16) 



If one has a resonance lying right at the threshold of 
channel 1 the couplings gi will be zero, since the deriva- 
tive of the loop function Gn, in the denominator of eq. 
(|16p is infinity at threshold. This is a general property 
which has its roots in basic Quantum Mechanics as shown 
in (2?| ■ In figure [1] we show plots of the real part of the 
loop function for the neutral and charged D meson chan- 
nels. 

It is interesting to note that eq. (JTHJ) for just one chan- 
nel is the method used to get couplings of bound states 




In the isospin symmetric case, the charged and neutral 
loops are equal, because these loops depend only on the 
masses, and therefore this ratio would be zero because 
the p contribution would vanish (no isospin violation). 

Actually the decays X — * J/ipp and X — ► J/ipu are 
not allowed because of phase-space, for p and lo with 
fixed masses, but can occur when their mass distribution 
is considered and will be seen in the decays X — > J/ipmr 
and X — > J/tJjtttttt respectively, where the two and three 
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pion states are the result of the decays of the p and u>. 
Hence to measure the ratio of the X decaying to two and 
three pions plus a J/ip one has to multiply the expression 



J 



in (|T7|) by the ratio of the phase-space available for the 
decay of a p to two pions divided by the phase-space for 
the decay of a a; into three pions: 



B(X -> J/i/jirir) 

B(X -> J/V-TTTTTr) 



Gn - G22 \ 2 Jo" V s ( s > m p, r p) g - rnj/^p - y/s) ds B p 
Gn + G22 J qS (s, m U) r w ) 6 (m x ~ mj/^ - </a) ds B u 

I 



(18) 



where B p and B p are the branching fractions of p decay- 
ing into two pions (~ 100 %) and u> decaying into three 
pions (~ 89 %), 9(y) is the Heaviside theta function and 
iS (s, m, F) is the spectral function of the mesons given 
by: 



S (s, m, r) = 



— Im 

IT 



1 



iTm 



(19) 



From the expression in eq. (|17p one observes that the 
isospin violation in the decay of the X will be propor- 
tional to the square of the difference between the loops 
with charged and neutral D mesons. Moreover if one 
looks at figure [1] one sees that this difference is maximal 
at the threshold of the D°D*°, such that the closer the 
resonance is to that threshold (the smaller the binding 
energy) the bigger is the isospin violation in the decay of 
the X. If the X is right over the threshold, the value of 
Rp/ U , with the loops calculated with dimensional regu- 
larization for p and lj fixed masses, is: 



R 



p/u 



0.032 



(20) 



This is a measure of the isospin violation in the decay 
of the X, which is only about 3% in spite of the fact that 
we have chosen the conditions to maximize it. This ratio 
is of the same order of magnitude as the one obtained in 
(see eq. (36) of this paper). However, even this small 
isospin breaking can lead to sizable values of the ratio of 
eq. (|18[) when one takes into account the mass distri- 
butions of the p and uj, which provide different effective 
phase-spaces in this two possible X decays. Thus, using 
eq. (|18p. which considers explicitly the p and u> mass 
distributions, we find the branching ratio: 



B(X -> J/iPn+n-TT ) 
B{X ^ JJ^t&tF) 



= 1.4 



(21) 



which is compatible with the value 1.0 ± 0.4 from exper- 
iment |6J|. 



III. THE TWO C-PARITY STATES 

There are six channels with charm and strangeness 
equal to zero and positive C-parity. We show in table 



TABLE I: Couplings of the pole at (3871.6-i0.001) MeV to 
the channels (a^=-1.27 here). 



Channel 


\gn-+pv\ [MeV] 


ir~ p + — c.c. 


1.4 


K~K*+ - c.c. 


8.7 


K°K*° - c.c. 


7.4 


D-D*+ - c.c. 


2982 


D°D*°-c.c. 


3005 


D~D* S + - c.c. 


2818 



[J the couplings of the pole obtained solving the scatter- 
ing equation for these channels. 

One can see in table 1 that, although there is some 
isospin violation in the couplings, it is very small, less 
than 1 %. One might think that if the binding energy 
is much smaller than the difference between the neutral 
and charged thresholds (8 MeV), the resonance will be 
mostly dominated by the neutral channel, the one clos- 
est to the threshold. The binding energy in the case of 
the pole in Table U is 0.4 MeV. As we mentioned, in the 
limit that the binding energy goes to zero, the couplings 
should all vanish. We show in figure [3] that, indeed, the 
coupling of the X to D°D*°goes to zero for small bind- 
ing energies, and in figure [4] we show that even though 
the difference between the neutral and charged couplings 
grows for small binding energies, they are of the same 
order of magnitude. The wave function of the A (3872) 
is, thus, very close to the isospin 1=0 combination of 
D°D*° — c.c. and D~ D* + — c.c. and has a sizable frac- 
tion of the D~D* + — c.c. state. 

From figure [4] we notice that indeed the isospin viola- 
tion in the couplings of the X to the DD* channels is big- 
ger for small binding energies, but it reaches a maximum 
of about 1.4% which is a very small value. We can go 
back to the argument that lead to eq. (flT]) and the only 
difference would be that the Gn and G22 functions would 
be multiplied by the D°D*° — c.c. and D + D*~ couplings 
from table |H which barely affect the results obtained in 
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FIG. 3: Coupling of the X to the D°D*° channel for different 
biding energies. 
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eq. (|21[) . since the differences in the couplings are much 
smaller than those between Gu and G 2 2- Although a 
D i> + c.c. state is proposed for the (3872) in p|, a 
formalism accounting for the ch arg ed component of this 
resonance is also presented in [38|. Our results would 
correspond to taking a value of the parameter 71 much 
bigger than «i(0) in size in eq. (39b) of '38]. However, 
no claims for any particular value of 71 are made in [3c| , 
where only the formalism is presented. 

As we already mentioned, we find a second state with 



negative C-parity. Some of the 22 channels with negative 
C-parity have isospin 1=0 to which the resonance can 
decay. There are also pure isospin 1=1 channels but, 
although the generated resonance is an isospin 1=0 state, 
these isospin 1=1 channels will couple to it since we are 
considering here some amount of isospin violation coming 
from the different masses of charged and neutral members 
of a same isospin multiplet. 

For values of an similar to those used in the genera- 
tion of the X(3872) (a.ff-1.27), the pole with negative C- 
parity is in the wrong Riemann sheet, but its effects can 
still be seen in the cross sections of some channels. We 
show in figure O the \T\ 2 plots of some channels. By tak- 
ing smaller values of an (around -1.36) one can recover 
a pole below threshold with Vs=(3871.4-i26.2) MeV. In 
our previous work of [10( this state was narrower. The 
reason for its relative big width in the present work is 
the inclusion of the 77-77/ mixing. As was explained in (34j 
the hidden charm dynamically generated states that we 
obtain are £(7(3) singlets and if one considers only the 
mathematical rj s without its mixing with a SU (3) singlet 
state 771 , the open channels for the resonance to decay are 
5(7(3) octets and are therefore suppressed. Only when 
considering also this singlet and hence the physical 77 and 
r)i states the open channels acquire a 5(7(3) singlet com- 
ponent to which the resonance strongly couples. 
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FIG. 5: The |T| 2 plot for some of the negative C-parity chan- 
nels. 



The channels shown in figure [5] are those where there 
is phase-space available for the resonance to decay and 
to which it couples most strongly. 

Next we want to compare the results obtained from 
our approach with the experiment data from [24| for 
B — > KD°D*°. For this we follow the approach of 
[l(il . [35j where one shows that the experimental data 
for dT/dM mv (DD*) are proportional to p\T DD , ^ DD ,\ 
where p is the center of mass momentum of the D meson 
with Mi nv invariant mass. 

We show in figure [5] plots of the dT/dMi nv for the 
channels D°D*° ± c.c. and the pure D°D*° and compare 
it with experimental data from (2~ij . 
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in spite of the results in the middle plot of figure [6l the 
experimental data does not rule out the existence of the 
negative C-parity state. 

The former discussions put the two states that we 
predict in a perspective concerning the D°D*° produc- 
tion experiment. In what follows we are going to do a 
more subtle exercise to bring some light into a current 
discussion on whether the combination of the data on 
X — > J/tfjirir and X —> D°D*° reactions determine if the 
state X(3872) is a bound state or a virtual one. In what 
follows we are going to consider only the contribution 
from the positive C-parity state. In [36| a slightly bound 
state is preferred, although a virtual state is not ruled 
out, while in [35[ a virtual state is claimed. With our de- 
tailed description of coupled channels, our approach is in 
a favorable position to get into the debate and bring new 
information. Yet, to do so one needs to introduce two 
new elements into consideration: the width of the D*° 
meson and the smearing of the results with experimental 
resolution. This was claimed to be relevant in [53] and 
[36| , We have considered this by taking for the D*° width 
T D ,a=65 KeV as in [36] and the experimental resolution 
AE=2.5 MeV. The consideration of the width of the D*° 
is taken into account by folding the D°D*° loop function 
G with the spectral function of the D* meson, 
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30 



20 
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-10 



n ft 313 
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12 



E [MeV] 

FIG. 6: The p|T| 2 , where p is the three momentum of the 
D mesons and E is the energy above threshold, compared to 
data from [3]. a = — 1.27 here. 



In the plots of figure [5] the theoretical curves have been 
normalized to fit the experimental data. One can clearly 
see that the positive C-parity state alone describes the 
data while the negative C-parity state alone does not 
describes it. However this experiment can not determine 
whether the D°D*° (together with the D°D*° summed 
incoherently) comes from a given C-parity In the lower 
plot of figure [5] we evaluate the differential cross section 
for the state D°D*°, which has contribution from both 
C-parity states. As we can see in the figure, the results 
obtained are also in agreement with the data and, hence, 



Im— — 

r: ) M 2 



Ml, 



iM D *T D * 



(22) 



in the calculation of the T-matrix, as done in eq. (20) 
of On the other hand the result for d ^ is folded 
with the mass distribution of the D* of eq. (|22"|) , since 
in the phase-space the three momentum q of the D°D* 
system appears as a factor and this three momentum 
depends on the mass of the D* . The final result is folded 
by a Gaussian distribution with a width of 2.5 MeV to 
simulate the experimental resolution. In this way one 
gets strength below the nominal threshold of D°D*° for 
the decay of the V(3872) into D°D*°. 

With these considerations we change slightly the a 
parameter which governs whether we obtain a bound 
state or a slightly unbound, virtual state. We normalize 
the two invariant mass distributions to the experimental 
data. The shapes alone tell us which option is preferable. 

In figures [7] and [5] we show the results for the dif- 
ferent values of a. To the left we have the results for 
X -> J/ipmr and to the right those for X -> D°D*°. 
What we see is that the effect of the convolution with 
the D* width and the experimental resolution is impor- 
tant, as claimed in [3(|[37j and help us make a choice of 
the preferred situation. At simple eye view, corroborated 
by a x 2 evaluation, see table [Til t ne preferred combined 
solution corresponds to a = —1.23 for which we have 
a slightly unbound, virtual state. This is the preferred 
solution in 135], also not ruled out in 36]. 
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TABLE II: \ 2 values for the fits of J/^n% and D° D*° produc- 
tion. The column to the right show the average value between 
the two. 





X 


X 




a 


for 


for 


x 2 






D°D*° 




-1.22 


1.83 


0.49 


1.16 


-1.23 


1.26 


0.50 


0.88 


-1.24 


0.87 


0.93 


0.90 


-1.25 


0.72 


2.77 


1.74 


-1.26 


0.92 


17.96 


9.44 


-1.27 


0.92 


20.31 


10.61 



IV. CONCLUSION AND OUTLOOK 

There is strong evidence that the AT (3872) state has 
J PC quantum numbers equal to It is very tempting 
to associate this state with a s-wave D°D°* molecular 
state if one takes into account the fact that the observed 
mass for this state is close to this threshold. 

Using a phenomenological Lagrangian and an unita- 
rization scheme for solving the scattering equation in 
couple channels, for all possible pairs of pseudoscalar and 
vector-mesons with zero charge and strangeness one ob- 
tains a pole with positive C-parity which can be associ- 
ated with the A (3872). Apart from this pole one also 
obtains strong attraction in the channels with negative 
C-parity, indicating the possible existence of a new state. 

Experimentally the decays of the A (3872) into J/tp 
with two and three pions have been measured to be of 
the same order of magnitude, suggesting a huge isospin 
violation in the decays of the A state. We have stud- 
ied here the effects of isospin violation in the decays of 
the A. The couplings of the A to charged and neutral 
D mesons are very similar, with at most 1.4% of isospin 
violation, in our model. As a consequence of that, once 
one considers the decay of A to J/tpp or J/ifsus pairs as 
going through DD* loops the J/ifsp production should be 
suppressed in relation to the J /ipus by about a factor 30. 
If one considers the decays of the A to J ftp with two and 
three pions as going through these channels respectively, 
one has also to take into account the phase-space avail- 
able for each decay which, due to the p width, is much 
bigger for this channel than for the us channel, compen- 
sating the factor 30 suppressing the decays of the A to 
J /t/sp- Moreover, the fact that the A is observed in these 
decay modes indicates a non negligible coupling of the X 
to vector-vector states, since it is visible in these channels 
(J '/isp and J/ipus) even with the small phase-space avail- 
able. This brought other authors to consider those chan- 
nels as possible extra building blocks of the A [3, ITsj ]. 
but we showed that we could explain the experimental 
phenomenology with dominant pseudoscalar-vector com- 
ponents. 

The prediction of a negative C-parity state is pecu- 
liar to DD* and would not appear as a vector-vector 
molecule. We predict this state with a framework which 



describes many low lying axial states and also most of the 
already observed axial charmed resonances. The com- 
parison of our theoretical invariant mass distribution for 
DD* production with data from Belle shows that one 
can not rule out the existence of this state and our model 
shows that the channels to which this resonance couples 
mostly are rj(f), T}us, r\ius and r) c us. We also made predic- 
tions for observables where the negative C-parity state 
should in principle be seen and we hope these results 
stimulate experimental efforts in this direction. 

Finally we made an investigation of the shapes of the 
A(3872) distributions in the J/ipirir and decays 
and found that they favor a slightly unbound, virtual 
state for this resonance. 
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APPENDIX A: THE £ COEFFICIENTS 
1. Charge Basis 

In charge basis we use the following ordering for the 
channels: 

1) p+71--, 2)K*+K-, 3)p°7T , 4)w7T°, 5)0tt°, 6) A, 
7)pV 8)K* a K°,_ 9)D*+D-, W)D* Q D°, ll)Ac, 
12)J/^7T°, 13)i^*°i^°, 14)p-7r+, 15V/, 16)0^7, 17)wt?/, 
18)077/, 19)L>* L> , 20)D*-D+, 21)ur) C) 22)07y c , 23)J/^r 1 , 
2i)J/ipV', 25)£>*+£>7, 26)DJ-D+, 27)K*-R+ and 
28)J/ipr) c 

Tables HIH HVl IV! PVTl and IVTTI show the coefficients £ 
for these channels. 



TABLE III: £ coefficients in charge basis basis for channels 1 
to 6. 



TABLE IV: £ coefficients in charge basis basis for channels 7 
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TABLE V: 
to 16. 
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£ coefficients in charge basis basis for channels 12 
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2. positive C-parity basis 



The channels here are the following: 
1)K* K° - c.c, 2)p+tt- - c.c, 3)D*°D° - c.c, 
4)D*+D- - c.c, 5)D*+D~ and 6)K*+R- - c.c. 

Table WTn\ shows the £ coefficients for these channels. 



3. Negative C-parity basis 



The channels here are the following: 

1) p+tt- + c.c, 2)K*_+K- + c.c, 3)p°n t) , 4)wtt°, 
5)07t°, 6)^°??, 7>V, 8)K*°K° + c.c, 9)D*+D- + c.c, 
10)D*°D° + c.c, ll)Ao 12 )J/i>K°, 13)wr), 14)077, 
15)wt?/, 16)0??/, 17)wt? c , 18)07/ c , 19)J/Vu 20)J/iM 
2l)D* + D- + c.c. and 22) J/^ 

Tables US El El and [XT]] show the £ coefficients for 
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TABLE VI: £ coefficients in charge basis basis for channels 17 
to 21. 



TABLE IX: £ coefficients in negative C-parity basis basis for 
channels 1 to 6. 
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TABLE VII: £ coefficients in charge basis basis for channels 
22 to 28. 
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TABLE VIII: £ coefficients in positive C-parity basis. 



i -> j 


1 


2 3 


4 


5 


6 


1 


2 


1 


7 


-7 


-1 


2 


1 


2 7 


7 





1 


3 





7 V + 1 


-1 


-1 


7 


4 


7 


7 -1 


4> + i 


1 





5 


-7 


-1 


l 


y> + 1 


7 


6 


-1 


1 7 





7 


2 



these channels. 



1 




2 


c 

? 










1 


2 


2 


2 


2 


1 


3 


i 


Q 
O 


9* /9 


A 
*± 


72 


A 
*± 


n 

u 


5 


-1 


5 





6 
7 

8 


2 


6 





v i 


7 



1 


V D 

1 


8 

Q 

£7 


9 





7 


10 


7 


10 


7 


11 





11 





12 





12 
13 






13 


2 
V 3 


14 





14 


-2^A 


15 


o 


15 


j/6 


16 





16 


73 


17 





17 





18 





18 





19 





19 





20 





20 





21 





21 


7 


22 





22 






3 




3 





4 





5 





6 





7 





8 
9 

1U 
11 


1 

72 

7 

72 
7 

73 




12 





13 





14 





15 





ID 


u 


17 





18 





19 





20 





21 





22 






4 


i 










-» 


5 


i 


6 


t 


4 





-» 




-> 




5 





5 





6 





6 





6 





7 





7 

8 
9 




1 


7 

8 




1 


8 
9 
10 


2 

V3 

7 


V2 

7 


9 

10 






7 

73 


10 


7 

72 


11 





11 


U 


i -i 
11 


U 


12 





1 9 


U 


1 9 
1Z 





13 





1 9 

lo 





1 9 

13 





14 





1 /I 

14 





14 


U 


15 





lo 


U 


15 





16 





16 





16 





17 





17 





17 





18 





18 





18 





19 





19 





19 





20 





20 





20 





21 





21 





21 





22 





22 





22 














TABLE X: £ coefficients in negative C-parity basis basis for 
channels 7 to 11. 
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TABLE XI: £ coefficients in negative C-parity basis basis for 
channels 12 to 16. 



TABLE XII: £ coefficients in negative C-parity basis basis for 
channels 17 to 22. 
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FIG. 7: Theoretical results compared to data from [2J|. The 
smeared points are calculated from the theoretical curve by 
folding it with a Gaussian, simulating the experimental reso- 
lution. 
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